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^ : Abstract 

Q-i! We study the power corrections (infrared renormalon contributions) to the coefficient 

Qh! functions for non-singlet deep inelastic structure functions due to gluon vacuum polar- 

ization insertions in one-loop graphs. Remarkably, for all the structure functions Fi, 
F2, -F3 and Qi, there are only two such contributions, corresponding to and l/Q^ 

power corrections. We compute their dependence on Bjorken x. The results could be 
used to model the dominant higher-twist contributions. 
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1 Introduction 



A number of recent papers have discussed power corrections to non-singlet deep inelastic structure 
functions and sum rules, using either the language of infrared renormalons or a dispersive 

approach based on assumed analyticity properties of the running coupling [^j. The contributions 
computed correspond to the sum of vacuum polarization insertions ('renormalon chains' [^) on the 
gluon line in one-loop corrections to the coefficient function. Phenomenologically, it appears 
that the corrections computed in this way provide a good guide to the form of the higher-twist 
contributions observed experimentally. This may be understandable in terms of the notion of a 
universal infrared- finite effective coupling [^]. 

In the present paper we point out that, within the framework of the approaches mentioned, 
there are only two power-behaved contributions, corresponding to 1/Q^ and corrections. This 
appears to be a special feature of deep inelastic structure functions. It is not the case for e~^e~ 
fragmentation functions, for example, where the same graphs generate an infinite set of l/Q^^ 
corrections. Of course, higher power corrections might also be generated by inserting renormalon 
chains in multi-loop graphs. Nevertheless, to the extent that the calculations based on one-loop 
graphs are phenomenologically successful, one may regard higher-power contributions as correction 
terms in the deep inelastic case. 

We present our calculation using the dispersive approach introduced in Ref. p. This makes 
it simple to compute the power-behaved contributions and their dependence on Bjorken x. As 
mentioned above, the results could be useful as a model for the dominant higher-twist contribu- 
tions. They are expressed in terms of two non-perturbative parameters, which could be determined 
experimentally by fitting deep inelastic data at moderate Q^. 



2 Dispersive method 

The method of Ref. starts from a (formal) dispersion relation for the QCD running coupling 
as(/c^) of the form 

a^ik^) = - psifi^) , p,(/i2) = -—Disc {a, (-/i^) I . (2.1) 

Introducing the effective coupling aes{fi'^), defined in terms of the 'spectral function' Ps{l^'^) by 

= Tr^«efr(/i') , (2.2) 
a m 

it follows that 

In the perturbative domain as ^ 1 the standard and effective couplings are approximately the 
same: 

«efr(/i') = «s(/i') - ^77"?^ + ... = «, + 0{al) . (2.4) 
o dm fi^ 



Thus one may regard defined by (|2.2| ) as an effective measure of QCD interaction strength, 
extending the physical perturbative coupling down to the non-perturbative domain. 
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Next we write the strong coupling in the form 



a,{k^) = a^''{k'')+Sasik^), (2.5) 

where a^"^ is the perturbative couphng and Sag is a modification to the effective interaction strength 
at small momentum scales, responsible for non-perturbative effects. The corresponding "effective 
coupling modification" da^s, which generates the non-perturbative interaction strength Sag, must 
satisfy a dispersion relation of the form 



S^sik^) = k'f ^^,%^, 5a..if^'). (2.6) 

It follows that an arbitrary finite modification of the effective coupling at low scales would generally 
introduce power corrections of the form l/k'^^ into the ultraviolet behaviour of the running coupling 
as itself. As discussed in Ref. P, such a modification would destroy the basis of the operator 
product expansion 0. One must therefore require that at least the first few integer moments of 
the coupling modification should vanish: 



/■CO djj^ / \p 

— (/i^j Sttesifl'^) = 0; p = 1, . . . ,Pmax • (2.7) 



J' 

The upper bound Pmax could be set by instanton-anti-instanton contributions at Pmax < /?o ~ 9. 

The effect on some observable F of gluon vacuum polarization insertions in one-loop graphs is 
represented in terms of the spectral function Ps(/i^) by a characteristic function jF(/i^), as follows: 

/i^ JO fi^ 



F = aMm + r ^Psip') ■ Hp') = r ^Psip') ■ [Hp') - m ] , (2. 

Jo u Jo u ^ J 



where we have made use of the formal relation (|2.1| ) to eliminate as{0)- The characteristic function is 
obtained by computing the relevant graphs with a non-zero gluon mass fi . Note that we do not 
intend to imply that the gluon has a real effective mass, but only that the dispersive representation 
( p.l| ) can be expressed in this way. 

Introducing the effective coupling a^sip'^) using Eq. ( p.2| ) and integrating by parts, we can write 

Fix,Q') = J^a,^{^^)-Hx,Q'\p'). (2.9) 

Here we have taken F to represent a quark structure function, with dependence on the deep inelastic 
scattering variables x and Q^. To obtain the corresponding (non-singlet) hadron structure function 
F(a;, Q^), we have to convolute the quark structure function with the appropriate combination of 
quark distribution functions q{x). As discussed in Ref. f^, one finds that the characteristic function 
JF has a coUinear divergent part, which generates the scale dependence of the quark distributions 
and the usual logarithmic scaling violations. The remaining part T'^'^^ generates the coefficient 
function C in the relation between the structure function and the quark distribution: 

/■I fl 2' 

Fix, Q') = / -Ciz, Q') qix/z, Q') (2.10) 

Jx Z 

where 

C{xM') = r%a,A^^)-P'\xM'\p')- (2.11) 
Jo 



Since T depends only on dimensionless ratios, we may write 

T{x,Q''-^^)=T{x,e) , T=-t^T{x,e), ^=^- (2-12) 

The power-behaved contributions depend on the small-e behaviour of J^, which is of the generic 
form 

J^{x,e) = -P{x)\ne + Coix)-C2ix)e\ne-lC4{x)eHne + --- . (2.13) 

The dots indicate terms that are analytic and vanishing at e = 0. Thus according to the above 
definition the regular part is 

r^^{x,e) = Co{x)-C2{x)e\ne-^C4{x)eHne + --- , (2.14) 

and 

J^''^{x,e) =C2{x)e\ne + C4{x)eHne + ■■■ . (2.15) 

The constraint (|2.7] ) means that only those terms in the small-e behaviour of the characteristic 
function that are non- analytic at e = will lead to power-behaved non-perturbative contributions 
0]. According to Eq. ( p.ll|) , the corresponding contributions to the coefficient function will be of 
the form 

C^P(a;, Q') = r ^ 6a,f,{f,')r'^{x, e = fi'/Q') . (2.16) 
Jo fi^ 

Thus from the small-e behaviour ( p.l5| ) we find 

C^P(x,Q^) = C2(x)^ + C4(x)^ , (2.17) 

where, following Ref. |p, we have defined the log- moment integrals 

A, = ^ ln(/.V/i^) 5aeff(/i^) . (2.18) 

Notice that since integer /i^-moments of ^ftefr vanish, these quantities are independent of the scale 
/Iq. For convenience, we extract a universal factor of Cpjln from the characteristic function. 

The interpretation of Eq. ( |2.17| ) in terms of a universal low-energy effective coupling is optional. 
More generally, we could interpret this expression as the ambiguity in the perturbative evaluation 
of the coefficient function, arising from the factorial divergence of the perturbation series generated 
by gluon vacuum polarization insertions. In this language, the two terms correspond to infrared 
renormalons, and the constants ^2 and A!^ are proportional to powers of the QCD scale A, the 
constants of proportionality depending on how we choose to resolve the renormalon ambiguity. 



3 Calculations 



In this section we apply the dispersive method to compute the power correction terms (|2.17 ) for 



non-singlet structure functions. Recall that the object of central importance is the characteristic 
function Tie) for the emission of a gluon with mass-squared /i^ = eQ^ at the hard scale . For 
power corrections, the relevant contributions are given by the non-analytic terms in the small-e 
behaviour of the logarithmic derivative ^(e). 



Q 



The characteristic function for the structure function F2 (actually F2/X) was given in Ref. 0: 
J^^{x, e) = J'f ^(x, e) 6(1 - X - xe) + J^^'"\t) 5{l - x) . (3.1) 



The contribution from real gluon emission is 



x, e 



2(1 -e)^ , 
1 -X ^ 

3 + 14e - 15e2 



1 + x) + 2(2 + X + 6x^)e - 2(1 + x + 9x^)e^ 
2 



In 



(1 — xe)(l — x) 



x^e 



e e 
+ 7^ ^ + 



+ 



X 



2(1 -x) (l-x)2 2(1 -x)3 l-xt 
l + 3x + 6(l-x)(l + 3x)e- 



+ 9x + 18x^)e^ . 



The virtual contribution is 

^('')(e) = 2(l-e)' 



where 



We note the relation 



Li2(e) + In e ln(l - e) - i In^ e - — 



Li2(e) = - r^ln(l-t). 
JO t 

^W(e) = - /V(")(x;e)e(l - X - xe)dx 

JO 



(3 - 2e) lne + 2e, 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



which means that the Adler sum rule is satisfied identically, that is, it receives neither perturbative 
nor power corrections (see Ref. |T0[|). 

Taking the small-x limit, we obtain an expression of the form (|2.13|) . The coefficient of — Ine 
in is the quark splitting function P(x) = (1 + x^)/(l — x), which is singular for x — > 1. The 
singularity is regularized by including the virtual contribution. As discussed above, this term pro- 
duces the usual logarithmic scaling violation. The second term Co(x) is the perturbative coefficient 
function (in the gluon mass regularization scheme). The remaining terms are analytic at e = 0, 
except for two terms proportional to elne and e^lne. Thus we obtain only two power corrections, 
of the form ( p.lTD . 

When taking the e ^ limit of eq. ( |3.2|) , we have to be careful with the functions that become 
singular in this limit. Defining '+', '++' and '+++' prescriptions such that, for any test function 
/, 



(3.6) 



F{x)+f{x)dx= F{x)[f{x)-f{l)]dx 
Jo 

^ F(x)++ fix) dx = /' F(x) [fix) - /(I) + (1 - x)f (1)] dx 
Jo 

F(x)+++ fix) dx = C Fix) [fix) - /(I) + (1 - x)f (1) - i(l - x) V"(l)] dx 



and recalling that 



5(")(1 -x)/(x) dx = /(")(1) 



(3.7) 



A 



we should replace the singular terms in Eq. ( p.2[ ) at small e, up to terms of order e^, as follows: 

+ e(l-e) 6'{l-x)-\e^ 6"{l-x) 



(e In e - In e - i In^ e - \e^) 5(1 - x) 
le2(lne-i)5"(l-x) 



1 — X \ 1 — X / 

N / -I- 

+ e(lne-elne- -x) 

^ 2 + '^(1 - ^) + e - e) - x) + ie(l - e) 

i X I L 



(1 -X)2 

fl -x)3 



■ X 



fl -x) 



+ (i + e) 5(1 - x) - e5'(l - x) - lne5"(l - x) 



Individual terms of the form e^lne with p > 2 are generated by the above replacements, but 
they cancel in the sum ( p.l|) . The only remaining terms in J^'^'^^ that are non-analytic as e — > are 
those of the form elne and Ine. 

The characteristic function for the structure function 2Fi is given by 

J='i{x,e) = J^2{x,e) - J^L{x,e) (3.9) 

where J-'l, the characteristic function for the longitudinal contribution Fl/x, is 



J^L(x,e) = 4(2 - 3xe)x^eln 



(1 — xe)(l — x) 



x^e 



2e(2 - e) ^ 



1 — X 



(l-x) 



(3.10) 



+ 2x + 4(1 - x)(l + 3x)e - 2(2 + 3x + 6x^)e^ . 
This contribution introduces no new power corrections, as noted in Ref. 0. 

The characteristic function for the parity- violating structure function was given in Ref. 0: 

^3(x,e) = JS(x, e) - :rrf(x,e) (3.11) 

where 



J^rf(x, e) = 2(4 + e - 9xe)x^e In 



(1 — xe)(l — x) 
x^e 



e(4 - 3e) 



1 — X 



(1 -x)- 



+ 



2x 



1 — xe 



(3.12) 



+ 1 - X + 2(2 + 5x - 9x^)6 - (5 + 7x + 18x^)6^ . 

Again, no power corrections beyond those in Eq. ( p.l7|) are introduced by this contribution. 

The power corrections to the coefficient function for the polarized structure function gi are the 
same as those for F^, since their characteristic functions are identical to one-loop order. 



4 Results 



The coefficients in Eq. (|2.171 ) for the power corrections to the coefficient function for F2/X are found 
from Eqs. (|3)-(|33) to be 

4 

C2(x) = — + 2(2 + X + 6x2) _g^n_^\_ _ 3,) 

(4.1) 

CAx) 



^ 4(1 + X + 9x3) ^ ^5 - x) + - 6"{1 - x) 

(l-x)+ 2 



The corresponding expressions in moment space, defined by 



C{N) = [\^~^C{x)dx (4.2) 
Jo 



are 

1,,2 3_ _ 4 4 36 



CJN) = -N' - -A^ + 16 4^1 



with 



N-l ^ 



For 2Fi = F2/X — Fl/x, the corresponding resuhs are[] 

C2(x) = - ^ + 2(2 + X + 2x^) - 5 5(1 - x) - 6\1 - x) 
(1 — xj+ 

CJx) = - — 4(1 + x + 3a;^) + 11 5(1 - x) + 4 5'(1 - x) + -5"(1 -x) 

(l-x)+ 2 



15 A A 12 

CAN) = -N^ + -N + 8-- ^ - ASi . 

^^'22 N N+l A^ + 3 



For F3 (and ^fi), 



4 

C2(x) = — + 2(2 + X + 2x2) - 5 (5(1 - x) - (5'(1 - x) 

(1 -x)+ 

C4(x) = ^ - 4(1 + X + x^) + 9 5(1 - x) + 4 5'(1 - x) + ^ 5"(1 - x) 
(1 — xj+ 2 

C2(N) = -iV - 4 + — + + + 4^1 

' N N + l N + 2 ^ 

-15 444 

CAN) = -N^ + -N + 6 4^1 . 

^^22 A^+1 A^ + 2 



(4.3) 



Si=T- = ij{N)+-fE = \nN + 0{l/N). (4.4) 
,=1 J 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



Note that the l/Q^ coefficients C2(x) are the same for F3 and 2Fi, but the l/Q^ coefficients C^i^x) 
are slightly different. 



5 Discussion 

To illustrate the above results we examine the l/Q^ and contributions arising from the valence 
quark distributions in neutrino scattering. Defining F = ^{F'^ + F'^)v and q = uy + dy, where V 
indicates the valence contribution, we can write 

F{x,Q ) ~ q{x,Q ) 1 + 



Q4 



(5.1) 

q{x,Q'^)Jx z 



A' /■! dz 
D2p{x, Q') = — ^ / -C2,{z) q{x/z, Q') 



^Our results for F]^ agree with those of Ref. 
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Figure 1: Coefficients of contributions to F2 (solid), and to Fi and F3 (dashed). 

the coefficient functions C2p being as given above for F = F2/X, 2Fi and F3. The valence quark 
distributions were taken from the MRSA parametrization [|1T|. Calculations were performed at 
= 10 GeV^, but the Q^-dependence of the coefficients D2p{x,Q^) is negligible. 

Figure shows the l/Q^ coefficients D2{x,Q^), assuming the value A'2 = —0.2 GeV^ for the 
non-perturbative parameter defined by Eq. (|2.18|) with p = 1. We see that this gives remarkably 
good agreement with the data points for F2, taken from Ref. |]12|-| The generally negative value 
of the prediction for F3 leads to a negative correction to the Gross-Llewellyn-Smith sum rule, in 
qualitative agreement with the predictions of Ref. 
CCFR-NuTeV Collaboration 



13], used in the test of the sum rule by the 



Figure ^ shows corresponding predictions for the 1/Q^ coefficients D^lx, Q^). Here the arbitrary 
choice = (^2)^ = 0.04 GeV^ is made for illustration only. We see that the contributions 
are peaked more sharply at high x, owing to the 6"{1 — x) term in C^lx). This indicates that the 
power corrections are functions of (1 — x)Q'^ rather than at high x. 

It would clearly be of interest to attempt a global fit to deep inelastic data at moderate Q^, 
treating the quantities A'2 and A'^ as free parameters. For this purpose, a program to compute the 
predicted power corrections is available from the authors. 
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